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Brady, Creighton and Thorne have proposed a choice of the lapse and shift for numerical evolu- 
tions in general relativity that extremizes a measure of the rate of change of the three-metric (BCT 
gauge). We investigate existence and uniqueness of this gauge, and comment on its use in numerical 
time evolutions. 



I. THE BCT GAUGE 



Cauchy data for general relativity consist of a three-metric g ab and extrinsic curvature K ab specified on a three- 
manifold ("slice") £. The Cauchy data determine the four-dimensional spacetime (locally) as a geometric object, but 
without fixing a coordinate system. When the spacetime is computed numerically as a sequence of spacelike slices 
£(i), the coordinates may be fixed incrementally by specifying a lapse a and shift vector (3 a on each new slice - a 
£f) • gauge choice. In one class of gauge choices, a and (3 a on a slice are determined by g ab and K ab on the same slice, i.e., 

o, 

O: (X,g ab ,K ab )^ (a,(3 a ). (1) 



Recently, Brady, Creighton and Thorne [JlJ have proposed a new gauge choice of this class in which a and (3 a are 
determined as a solution of the coupled equations 

^ K ah F ab = 0, D a F ab = 0, (2) 

, where D a is the covariant derivative with respect to the three-metric g abl and 

•i-H . 

X ! F ab {a, 13) = L ab (f3) - 2aK ab , L ab (/3) = D a /3 b + D b f3 a (3) 



These two equations arise when one varies the action principle 

1 = J F ab F ab ^d 3 x (4) 

with respect to a and /3 b , respectively. In a four-dimensional context, F ab is the time derivative of the three-metric 
induced on the time slicing with lapse a and shift (3 a : 

g ab = C t g ab = F ab (a, t a = an a + /3 a , (5) 

where n a is the unit normal vector of the slice £ when it is embedded into a spacetime. Equations (^|) are therefore 
also called the minimal strain lapse equation and minimal strain shift equation, where / is the "strain" that is being 
extremized. 

The motivation for considering these equations is that a good gauge choice should have the property (among 
others) to be compatible with approximate Killing vectors, in the sense that if an approximate Killing vector exists, 
the spacetime metric in that gauge should evolve as slowly as possible. The inspiral phase of a binary black hole 
system in a spherical orbit, for example, has an approximate Killing vector, and one would like to be able to find 
(corotating) coordinates in which the spacetime metric evolves on the timescale in which the orbit shrinks through 
the emission of gravitational waves, rather than on the much shorter orbital period timescale. 
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The first of the two equations (g) can be solved algebraically for a, 

K cd L cd (P) 



K e fK, 



ef 



(0) 



and the result substituted into the second equation. One obtains a vector linear second-order differential equation for 
the vector (3 alone: 

D a H ab ((3) = 0, H ab ((3) = L ab {(3) - 2^^K ab . (7) 

We shall call this single vector equation the BCT equation, and in the following we shall consider the lapse a as a 
dependent quantity determined by (^3j). Equation (^) can be obtained from the action principle 

J = J H ab H ab ^d 3 x, (8) 

which is obtained by substituting the lapse (Q) into the action principle /. 



II. EXISTENCE AND UNIQUENESS 



If one is to use the BCT gauge, it is important to know the answer to the following question: For which choices 
of data (U, g ab , K ab ) on a region U of £ with boundary dU and for which sets of boundary conditions for (3 can one 
solve the BCT equation (]?])? In this brief note, we show that for some choices of data there is a unique solution with 
specified Dirichlct boundary conditions, while for others, there are many solutions with those boundary conditions. 

The issues of existence and uniqueness of the BCT gauge have been previously considered by Gongalves [Q. 
Gongalves shows that the differential operator defined by (Q) is strongly elliptic if and only if K a b , considered as 
a map on the tangent space of E, has at most one vanishing eigenvalue. (Let the principal part of the differential 
operator acting on (3 a be M a b cd D c Ddf3 b . The operator is then defined to be strongly elliptic if M a b cd is positive 
definite with respect both to the two indices that slot into derivatives and the index pair that shuffles the vector index 
on (3.) Generally, operator ellipticity - strong or otherwise is not enough to determine whether a boundary value 
problem admits a solution. However, using the fact that the equation of interest (Q) is of divergence form, Gongalves 
does obtain local results. In this paper, using the Fredholm alternative, we obtain stronger, global, results. 

Before deriving the condition for existence of the BCT gauge, we wish to turn the boundary value problem for 
Equation (^) into one with homogeneous boundary conditions. Let us write Equation (j?]) in operator form as O(0) = 0, 
and the corresponding Dirichlet problem as 

0(0) = 0, P\ BU = f (9) 

for some given continuous vector-valued function /. Presuming that the region U and its boundary are well-behaved, 
we may always extend / to a function F : U — > i? 3 which is C 2 on U and continuous on U U dU . Then if we find a 
solution £ of the boundary value problem 

O(0 = -O(F), $\au = 0, (10) 

and if we set (3 — £ + F, we have f3 satisfying the boundary value problem (^). We may now focus on the discussion 
of the boundary value problem ( [To| ) (for arbitrary F) . 

To be able to use Fredholm ideas to study the boundary value problem (|Io|), one first needs to establish that (|Io| ) 
defines a Fredholm map. As shown in propositions 11.10 and 11.16 of ||, this fact follows from the strong ellipticity 
of O. We also need the following: 

Lemma: For (3 satisfying the Dirichlct condition f3\gu — 0, O is a self-adjoint operator. 

Proof: We consider the quantity 

/ iO{fl)= [ j a D b H ab (P) (11) 
Ju Ju 

with H ab from Equation (Q). Integrating by parts, and using ~i\ou = 0, we obtain 
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1 0{(i) = - / DP-fHM = - / F Q0 ( 7 )-Da/3 6 . (12) 
u Ju Ju 



Integrating by parts again, and using (3\ou = 0, we have 



7 0(/3)= / D a H ab ( 7 )p b = / /30( 7 ). (13) 

Hence O is self-adjoint on functions /3 that vanish on the boundary. (As the boundary term in the integration by 
parts is of the form J gu s /3 f, i? a h(7), where s° is the unit normal vector to dU, the operator O is also self-adjoint on 
the space of functions 7 that satisfy s a i/ a fc(7) — on the boundary. This boundary condition is the mathematical 
analog of the Neumann boundary condition for the Laplace equation, but it contains additional terms that make its 
physical meaning unclear. Therefore we do not consider it here.) 

Since O is self-adjoint, and since it defines a Fredholm map, the Fredholm Alternative ||] specifies a clear condition 
for existence of solutions. Stated for the operator O, one has 

Proposition (Fredholm Alternative): Fix (U, g ab , K ab ). The Dirichlet boundary value problem jl0| ) has a unique 
solution for every choice of O(F) if and only if the only vector function satisfying 

0(0-0, Z\ BU = (14) 

is £ a = 0. (That is, the kernel of O is empty.) For a given O(F), the boundary value problem (|l^) has a solution so 
long as O(F) satisfies 



f £O(F)=0 (15) 
Ju 



for every £ which satisfies (Q). (The kernel of O is orthogonal to the source.) 

We base our existence results for the BCT gauge on the Fredholm Alternative. So we are led to consider solutions 
£, of (|lj), or elements of the kernel of O. Given such a solution, we consider the quantity 

J J 0(0 = -~ J^iOHabiO = -\ ^H^iOHabiO, (16) 

where the first equality follows from integration by parts (£ vanishes on dU) and the second equality follows from 
K ab H ab = 0. £ therefore satisfies (0) if and only if 

Hab{0 = 0, £\ du = 0. (17) 

Thus the kernel of O consists of solutions of this boundary value problem. 

For a given set of initial data (U, g a b, Kab), the Fredholm Alternative asks that we determine the kernel of O. If the 
kernel is empty, then the Dirichlet boundary value problem (|Io|), and hence (0), admits a unique solution. 

If the kernel is not empty, with non-trivial elements £, then to determine the solubility of the boundary value 
problem ( |l0|) we need to consider J u £ i O(F). Integrating by parts, we have 

JjO{F) = -\ J^L{i)H{F) = -\ jf ff(0L(F), (18) 

but H(£) = by assumption. So this always vanishes, and a solution always exists in this case, too. It is determined, 
however, only up to the addition of any element of the kernel. 

We conclude that a solution of the BCT gauge exists whenever O is strongly elliptic, but if the kernel ( |l7j ) is not 
empty, the solution is not unique. 

III. EXAMPLES FOR UNIQUENESS AND NON-UNIQUENESS 

It is useful to note that there are choices of data (U, g a b, K a b) for which each of these two cases hold. For the first 
case, where the kernel is empty, we consider data with K a t = pg a b, with p nowhere vanishing. The kernel equation 
( |l7| ) then becomes 

Datb + D b £ a - ^gabD c £ c = 0, (19) 
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This is the equation for a conformal Killing vector. If we choose a metric g a t, that does not admit a conformal Killing 
vector that vanishes on the boundary, we have constructed data that give rise to an empty kernel, and therefore a 
unique BCT gauge for a given choice of the boundary data. 

For the other case, we consider a slice through a spacetime that has a timelike Killing vector. We choose the lapse 
and shift so that t a is the Killing vector. The time derivative of the three-metric then vanishes, and consequently 
F a b{ct,(3) — 0. From this it follows that H a b{(3) = 0. In order to obtain (3\qu — 0, we choose the slice so that it is 
normal to the Killing vector in dB but not in B. If in such a situation one tries to solve the BCT equation with a 
boundary where the slice is approximately normal to the Killing vector, the numerical problem might become badly 
conditioned. 

As a concrete example, we consider a spherically symmetric slice in flat spacetime. Let (t, r, 6, ip) be the standard 
spherically coordinates on Minkowski spacetime, and let the slice be given by t — T(r). As coordinates intrinsic to 
the slice we use (r, 9, if) induced by the spacetime coordinates of the same name. The normal vector to the slice, 
induced 3-metric and extrinsic curvature are then given by 

rpl -, 

(20) 



VI - T' 2 ' VI - T' 2 ' 

g„. = 1-T' 2 , g e0 = r 2 , g w = r 2 sin 2 6, (21) 

T" rT' 
Krr = ~^/T=W Kee ^~^T^W K vv = S m 2 eK B e. (22) 

The normal vector is parallel to the Killing vector d/dt where T'(r) = 0, say at r = r$. We choose U to be the ball 
r < tq. The desired element of the kernel of O is then 



a = Vl -T' 2 , (3 r = —T', (23) 

up to an overall constant factor. 

The potential difficulty with non-uniqueness can be avoided by an appropriate choice of slice and boundary. In 
particular, for the choice of slice and boundary proposed in H for the binary black hole problem, the shift is nowhere 
small, so that the Killing vector is nowhere normal to the slice. Near the black hole excision boundary, the slicing is of 
the Painleve-Gullstrand (or Kerr-Schild) type, with a large radial shift, while at the outer boundary the coordinates 
are corotating, with a large d/d<f> shift component. 



IV. USE OF THE BCT GAUGE AND OTHER GAUGES IN TIME EVOLUTION 



Consider the evolution of initial data (g a b, Kab)- The evolution equation for K a b is of the form 



K ab 



-D a Dba + other terms. 



(24) 



The time derivative of K a i, contains the second spatial derivative of the lapse. If one uses the BCT gauge to determine 
a and /3 from g a b and K a },, the evolution equation for K a f, becomes 



K nh = 



L cd {(3) 

K ef Ref 



l (P)K mn K 



cd 



(K ef K*f) 2 



D a DbK c d + other terms. 



(25) 



This means that if K a {> is initially in a function space of finite differentiability, time evolution takes it out of that 
space - it "loses two derivatives" . ((3 a itself appears to gain two derivatives because it is the solution of an elliptic 
equation, but this does not affect the argument.) 

While this is a technical obstacle to proving existence and uniqueness of solutions to the Einstein equations in the 
BCT gauge, it also hints at the possible existence of a practical problem for the use of the BCT gauge in numerical 
evolution. Roughly speaking, one would expect numerical noise to be amplified during time evolution, whereas 
maximal slicing with zero shift, for example, is empirically known to dampen numerical noise. In the toy model 
equation it = ku xx , noise is damped for k > (heat equation) and increases for k < (heat equation run backwards), 
so that one only needs to choose the correct sign of k. In contrast, Equation ( p5| ) is a nonlinear tensor equation. It 
appears plausible that some of the eigenvalues of its linearization around some backgrounds (g a b, K a b) correspond to 
negative K, thus leading to the growth of linear noise. Corresponding nonlinear instabilities may also exist. We have 
not attempted to investigate this question. 
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Here we would like to draw attention to an alternative gauge choice. Consider the following coupled equations for 
the lapse and shift: 

D b D b (i a + D b D a (3 b - 2D b (aK ab ) - ~D a (2D b l3 b - 2aK) = 0, (26) 



D a D a a 



a + (3 a D a K = 0. (27) 



Here r and p are matter terms, namely the trace of the (three-dimensional) stress tensor and the energy density, n is 
a constant that is either zero or one. These equations for a and (3 are elliptic, but they are not self-adjoint. Existence 
and uniqueness of solutions will be considered elsewhere As the equations are elliptic in both a and /?, (g ab , K ab ) 
do not lose differentiability compared to the Cauchy data (g a b, K ab ). 

The first, vector, equation is D b F ab = for n = 0, and is D b (F ab - l/3g ab F c c ) = for n = 1. These two shift 
conditions were suggested by Smarr and York in their classic paper on coordinate conditions S under the names 
"minimal strain shift" and "minimal distortion shift". The second, scalar, equation is K = 0. For K = it reduces 
to maximal slicing. Here, however, we do not assume that K is zero, nor that it is constant in space. K = slicing, 
and its combination with minimal distortion shift, was suggested by Smarr and York, and they also discuss its Killing 
vector-tracking property. Their subsequent discussion focuses on maximal {K = 0) slicing with minimal distortion 
shift. This more restricted gauge choice is now commonly associated with the name "Smarr-York" (SY) gauge, while 
the more generalgauge with K ^ seems to have been forgotten. The desirability of Killing vector tracking was later 
rediscovered in jij and || , and the gauge discussed here was rediscovered as "generalized Smarr-York" (GSY) gauge 
in§. 

Here we want to point out that the GSY gauge has the desirable properties of the BCT gauge - it tracks Killing 
vectors and it admits generic Cauchy data - without having the "loss of derivatives" property and the possible problem 
when the slice is normal to the Killing vector at the boundary. In a direct numerical comparison of the BCT and GSY 
gauges in spherical symmetry, Garfinkle and Gundlach || find that the GSY gauge stably tracks Killing vectors, but 
are unable to obtain a stable time evolution with the BCT gauge, and this may be due to the "loss of derivatives" 
property. Shibata reports good experiences with an approximate implementation of maximal slicing with minimal 
strain shift in 3D evolutions of a neutron star binary. Independent numerical work, and tests in 3D, are required to 
decide if one of the gauges is more suitable in practice than the other. 
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